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ABSTRACT 


Most  commercial  and  military  systems  subject  to  failure  have 
an  initial  period  of  growth  With  uncertain  time  dependent  growth 
rate.  Prediction  of  spare  parts  demand  and  service  personnel  demand 
is  correspondingly  uncertain  and  existant  statistical  tools  are 
inadequate  for  the  adaptive  ad  hoc  planning  needed. 

In  our  model,  systems  subject  to  failure  enter  into  use  at 
the  epochs  of  a  time-inhomogeneous  Poisson  process  of  rate  X(t). 

A  component  or  module  of  each  system  in  use  has  constant  failure 
rate  p  and  generates  demand  for  parts  and  service.  The  distribution 
of  the  cumulative  failures  N(t)  is  obtained.  Numerical  methods 
and  the  asymptotic  distribution  for  large  t  are  described. 
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In  Section  2,  the  case  of  constant  Poisson  build  up  where  X(t)  «  X 
is  treated.  The  cumulative  demand  distribution  at  time  t  is  obtained 
through  generating  function  methods  by  treating  the  population  of  systems 
in  use  and  the  cumulative  demand  as  a  bivariate  Markov  chain.  An 
algorithm  for  finding  the  distribution  numerically  when  the  cumulative 
demand  at  time  t  is  modest  is  displayed  in  Section  3.  A  graphical 
display  shows  growing  normality  for  large  demand.  In  Section  4,  time 
dependent  entry  rate  X(t)  is  treated,  and  the  distribution  of  cumula¬ 
tive  demand  N(t)  is  found.  It  is  shown  that  this  demand  is  both  a 
compound  Poisson  distribution  and  a  mixture  of  Poisson  distributions, 
and  that  the  mean  and  variance  of  the  demand  at  time  t  have  simple  useful 
forms.  Specifically  it  is  shown  that 

t 

E[N(t)]  =  W  /  X(T)(t-T)dT  (1.1) 

0 

and  that 

,  t 

Var[N(t)]  =  li  /  X(T)(t-T)^dT  ♦  E[N(t)].  (1.2) 

0 

The  asymptotic  character  of  the  distribution  of  demand  is  exhibited  in 
Section  5,  and  conditions  under  which  the  cumulative  demand  is  normal 
or  Poisson  are  described.  Practical  use  of  the  tools  developed  in  this 
paper  is  discussed  in  the  final  section. 

The  paper  is  self-contained.  Only  a  knowledge  of  probability  theory, 
Poisson  processes,  the  central  limit  theorem  and  elementary  stochastic 
process  tools  as  found,  e.g.,  in  Ross  (1970)  is  required. 
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!.  Replacement  Demand  Associated  with  Constant  Poisson  Build  Up 

Consider  a  random  counting  process  M(t)  for  "primary"  events  (e.g., 
aircraft  entering  service).  We  suppose  throughout  this  section  that  all 
primary  units  produced  are  still  functioning  at  time  t.  Each  primary  event 
generates  subsequent  secondary  Poisson  events  (e.g,,  demands  for  service 
or  part  replacement)  with  rate  y.  If  T.  is  the  time  of  the  i-th  primary 
event  and  if  L(t  -  T^)  is  the  random  number  of  secondary  events  up  to  time 
t,  initiated  by  the  primary  event  at  ,  then  the  total  number  N(t)  of 
events  is 


Mft) 

N(t)  =7  L(t  -  T.)  .  (2.1) 

i  =  l  ^ 

with  the  usual  convention  that  an  empty  sum  has  the  value  0.  We  are  inter¬ 
ested  in  the  distribution  of  N(t)  when  M{0)  *  0. 

We  first  consider  the  case  for  which  the  primary  counting  process 
M(t)  is  Poisson  with  rate  X.  To  find  the  distribution  of  N(t)  for  this 
case  one  can  use  bivariate  Markov  processes  and  bivariate  generating  func¬ 
tions.  Let  [M(t),  N(t))  be  the  bivariate  Markov  process  on  the  set  of 
states  {(m,n);  0  s  n,  0  s  n)  with  transition  rates 

'’(m,n)  (m*l  ,n)  “  ^ 

and  (2.2) 

'’(m,n)(m,n+l)  *  ’ 

Let  „(t)  -  P[M(t)  ■  m,  N(t)  ■  n|M(0)  ■  0,  N(0)  ■  0]  and  let  us  denote  by 

in  I  Ti 

g(u,v,t)  ■  I  p_  _(t)u"'v"  (2.3) 

all  m,n 
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the  bivariate  probability  generating  function  of  [M(t),  N(t)]. 
Since 


d  p  (t) 
^m,n  ’ 

dt 


(2.4) 


one  finds 


=  -X(l-u)g(u.v,t)  -  liu  1^  g(u.v.t) 

'  ““V  1^  g(u,v,t)  (2.5) 

If  g(u,v,t)  =  exp{Q(u,v,t) }  then  clearly,  for  g(u,v,0)  *  1,  one  has 

^CjCUjV.t)  ^  .x(i.u)  -  uu(l-v)|^  Q(u,v,t)  .  (2.6a) 

Q(u,v,0)  =  0  (2.6b) 

If  one  seeks  a  solution  linear  in  u,  i.e., 

Q(u,v,t)  =  A(t)  ♦  B(v,t)u  (2.7) 

then  (2.6)  becomes 

^  A(t)  4  u  .  .x(i.u)  -  u(l-v)u  B(v,t)  .  (2.8) 

One  then  must  have 

3j*<t)  .-X 


1^  BCv.t)  »  A  -  w{l.v)B(v,t)  . 


One  thus  obtains  a  particular  solution  of  (2.6)  with  A(0)  ■  B(0)  =  0: 


A(t)  *  -At  , 


and 


B(t)  =  Xt(i 


g-yt(l-v) 

lit  (1-v)  ^ 


Hence 

g(u,v,t)  =  exp{-Xt[l  -  .  (2.9) 

A  simple  argument  demonstrates  that  (2.6)  has  only  one  solution,  which  we 
have  found. 

The  probability  generating  function  for  N(t)  is  then 

o^(v)  »  =  gd.v.t)  (2.10) 

which  becomes 


-vt(l-v) 

P^(v)  =  exp{-Xt[l - ut'(l-vj - •  (2.11) 

This  is  the  generating  function  of  a  compound  Poisson  process,  whose 
increments,  say  V^,  are  simply  related  to  a  secondary  Poisson  distribution. 
Specifically,  one  can  write  (2.11)  as 

c^(v)  =  exp(-Xt[l  -  o^(v)]>  (2.12) 


where 
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(v) 


E[v  = 


,  _  -Ut(l-v) 
Pt(l-v) 


(2.13) 


o^(v)  is  then  the  generating  function  of  the  compounding  distribution  of 
V^.  It  is  easily  seen  that  a^(v)  corresponds  to  the  distribution 


P(V^ 


1 

I 

r=n*l 


(wt)^  -ut 
r!  ® 


(2.14) 


The  mean  and  the  variance  may  be  obtained  directly  from  (2.11)  and  are 


E(N(t)] 


(2.15) 


Var[N(t)] 


(2.16) 


Note  that  the  variance-mean  ratio  is 


VarfN(t)1 

EtN(t)] 


(2.17) 


which  is  always  larger  than  one. 

We  note  that  the  generating  function  (2.11)  is  similar  in  character 
to  that  of  a  generating  function  obtained  by  Berg  (1981). 
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>.  Evaluation  of  the  Distribution 

We  have  just  shown  that  the  distribution  of  N(t)  is  that  of  a  com¬ 
pound  Poisson  distribution  with  the  underlying  Poisson  process  having 
rate  X,  and  with  the  compounding  variable  having  the  distribution  (2.14). 
Thus, 


P_(t)  =  P{N(t)  =  n)  =  J  e*^'  a^t)  (3.1) 

k=0  n 

k 

where  “j^(t)  is  the  n-th  term  of  the  k-fold  convolution  of  with  itself. 
In  other  words, 

a^(t)  =  P{V*’^  =  n)  .  (3.2) 

*k 

where  denotes  the  sum  of  k  independent  random  variables,  each  having 
distribution  (2.14).  The  distribution  of  can  be  obtained  easily  by 
convolving  iteratively  the  distribution  (2.14). 

One  can  obtain  Pj^(t)  explicitly  from  (2.11)  for  small  values  of  n 
by  differentiation.  One  finds  that 


Po(t) 

PjCt) 


exp(-Xt  ^  (1  -  e'‘'^)) 


3o^(v)| 


9v 
-  wt  e 


v»0 

vt 


1 


(3.3) 


(3.4) 


and  that 


^2^'^  *  J  P.(v) 
3v^  ^ 


v»0 


♦  2(1  -  e'^')  -  2ute‘‘''  -  (ut)V''*} 


(3.5) 
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For  larger  values  of  n,  say  n  i  3,  the  probabilities  P^jCt)  can  be  obtained 
numerically.  The  plots  of  the  standardized  survival  functions,  i.e.,  of 

p^N(t]  -  EfN(t)l  ^ 

Aar[N(t)] 

for  various  values  of  parameters  B  =  Xt  and  A  =  pt  on  normal  probability 

r 

paper  are  shown  in  Figure  1.  The  diagonal  straight  line  in  Figure  1  repre¬ 
sents  the  survival  function  of  the  standard  normal  distribution. 

From  the  plots,  it  can  be  seen  that,  for  about  A  i  20  and  B  i  10, 
the  central  part  of  the  distribution  can  be  well  approximated  by  the 
standard  normal.  Note,  that  the  left  tail  of  the  distribution  of  N(t) 
is  heavier  than  normal  so  that  the  normal  approximati.>n  underestimates 
the  percentage  points  of  the  N(t)  distribution.  On  the  other  hand,  the 
right  side  of  N(t)  is  lighter  and  the  mass  point  at  the  extreme  left 
represents  the  probability  of  having  no  replacements  needed.  Theoret¬ 
ical  justification  of  the  normality  will  be  given  in  Section  5. 


mil 
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The  Distribution  of  the  Cumulative  Demand  NCt)  for  More  General 
Primary  Input  Streams 

The  process  described  in  Section  2  can  be  generalized  to  better  model 
real  world  situations. 

Theorem  4.1.  If  the  sequence  of  primary  events  is  any  point  process  and 
if  the  number  of  secondary  events  generated  by  primary  events  is  Poisson, 
then  NCt)  is  a  mixture  of  Poisson  distributions,  i.e., 

P*(v)  =  =  /  e‘®^^'''^dG  Ce)  (4.1) 

0 

where  0^(6)  is  a  c.d.f. 

Proof :  Let  oi  =  (t ^ (u) (w) , . . . )  be  a  sample  sequence  of  primary  epochs. 
Then  the  conditional  probability  generating  function  is 

-a^  (u) (1-v) 

and  a^(ui)  is  a  Poisson  parameter.  Hence 
P^(v)  =  /  P^(vjja^(dw)  , 

where  o^Cdu)  is  a  measure  on  the  space  of  sample  sequences.  The  theorem 
then  follows.  0 

Example  4.2.  We  will  verify  the  Poisson  mixture  property  for  the  Poisson 
input  case  of  Section  2.  For  notational  brevity,  let  it  =  B  and  ut  =  A. 
When  the  primary  stream  is  Poisson,  (2.11)  may  be  written  as 


p^  (v 


-11- 


Letting  1-v  ^  s,  ve  get 

/  e  ”^dG(e)  =  expf-B[l  -  1.:-® - ]}  . 

0 

But  this  is  a  Laplace  transform  of  a  compound  Poisson  process  obtained 
from  a  uniform  distribution  over  (0,A).  Hence  the  mixing  distribution 
G(6)  is  that  compound  Poisson  distribution. 

Remark  4.3.  The  result  of  Theorem  4.1  is  not  altered  by  phase  out  of 
primary  sources  due  to  failure,  obsolescence,  or  random  utilizations, 
as  we  see  next. 

For  most  industrial  applications,  the  primary  production  will  also 
vary  with  time  due  to  initial  buildup  and  ultimate  phaseout.  Not  only 
will  primary  units  (e.g.,  aircraft)  disappear  because  of  catastrophic 
failure  or  obsolescence  but  the  primary  units  available  may  have  random 
utilization  arising  from  economic  fluctuations,  long  downtimes,  etc. 

For  such  a  more  general  input  streams  let  M(t)  be  the  total  number 
of  primary  units  in  operation  at  time  t,  and  let  N(t)  be  the  cumulative 
number  of  replacement  parts  (service  procedures)  required  up  to  time  t. 
Let  m^(t)  be  a  realization  of  the  input  process.  For  each  such  realiza¬ 
tion  let  N^(t)  be  the  random  variable  associated  with  that  primary 
realization.  If  the  secondary  events  are  Poisson  with  rate  y,  then 
N^(t)  has  Poisson  distribution  with  parameter 

t 

Iy(t)  -  V  I  my(y)dy  .  (4.2) 

Equation  (4.2)  defines  a  random  process  I(t)  of  basic  interest.  Since 
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I(t)  *  (wt)  •  i  /  M(y)dy 
0 


(4.2a) 


we  see  that  I(t)  is  dimensionless  and  that  (wt)‘^l(t)  is  the  average 
population  up  to  time  t.  Moreover, 


-I  (t)Cl-v) 

P^(v)  =  E[v  ^  ^]  =  /  e  a^(dui) 

.E[e-^(t)(l-v)3  ^ 


(4.3) 


Comparing  (4.1)  and  (4.3)  we  see  that 


G^(e)  =  P{I(t)  <  e)  . 


(4.4) 


The  following  results  are  a  direct  consequence  of  Theorem  4,1. 


Theorem  4.4.  (a)  E[N(t)l  *  EfKt)!  =  w  /  E[M(y)]dy 

0 

(b)  E[N(t)(N(t)-l))  =  E[I^(t)] 


(4.5) 

(4.6) 


(c)  Var[N(t)]  *  VartKt)]  ♦  E[I(t)]  .  (4.7) 

Note  that  if  I(t)  is  non-random  with  E[I(t)]  =  E(t),  then  we  have 


(4.8) 


and  thus  N(t)  is  a  nonhomogeneous  Poisson  process  with  rate  E(t),  i.e.. 


p  (t)  .  ^"C(t)  ^  ,  0,1,2,.., 


(4.9) 


Example  4.5.  Suppose  that  the  primary  process  is  deterministic  and  the 
interarrival  times  between  primary  events  equal  T  .  In  this  case 


ni^(y)  ^  k  ,  if  kT*  s  y  <  (k^l)T* 
for  k  *  0,1,2,...,  and 


where  K  =  [t/T*]. 

Theorem  4.6.  If  the  primary  production  stream  is  a  nonhomogeneous  Poisson 
process  M(t)  with  rate  X(t)  continuous  in  t  then  the  number  of  secondary 
events  N(t)  has  a  compound  Poisson  distribution. 

Proof :  We  have  seen  in  (4.2)  that  N(t)  is  a  Poisson  mixture  with 

random  mixing  parameter 


t 

l(t)  «  g  /  M(T)dT  =  gS(t)  ,  (4,11) 

0 

t 

where  S(t)  »  {  M(T)dt  is  a  stochastic  integral.  One  then  has  to  first  order 
0 

in  dt, 

Efe*sdM(T)j  ^  ^-X(T)dtll-e'*] 
and 


where 


♦s(t)''>  ■  El'”®"’) 


/ 


I 
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t  y  t 

Hence,  S(t)  •  [  dy  J  dM(T)  ■  /  (t-T)dM(T)  for  every  sample  path.  For 
0  0  0 

any  Riemann  sequence  of  partitions  of  the  interval  [0,t)  whose 

maximiim  interval  6.  *  max  6..  goes  to  zero 
3  k 

S(t)  -  lim  I  (t-T.)(M(T  )  - 

6^-»0  k 


-(t-T.  )s 

”  exp{-X(T.)  6..  [1-e  ]} 


-  lim  exp{-y  X(t.)  «^v[l-e  ]) 

6  -►0  k  »•  jx 
3 

*  exp{-/  X(T) (i-e‘ 

O 


«  exp{-A(t)tl  -  /  ^  dTl) 


where  A(t)  ■  /  X(T)dT.  Thus  (4.12)  becomes 
0 


♦^(v)  ■  exp{-A(t)[l  -  /  di]} 


(4.13) 


exp{-A(t)(l-a^ (v)]}. 


(4.14) 


This  is  the  generating  function  of  a  compound  Poisson  distribution,  with 


a.,v,  .  I  .-(.-Ov(l-v)  ^  „  , 


(4.15) 


Remark  4.7.  For  the  homogeneous  Poisson  case  with  X(t)  ■  X,  (4.15) 
reduces  to  a^(v)  of  (2.13). 
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5.  Asymptotic  Normality  and  Asymptotic  Poisson  Character  of  the 
Cumulative  Demand 

Consider  first  the  time  homogeneous  process  of  Section  2.  When 
the  parameters  A  «  pt  and  B  «  Xt  are  small,  one  expects  that  the  distri¬ 
bution  of  N(t)  will  be  Poisson.  One  also  expects  that  when  A  and  B 
are  large  the  distribution  of  N(t)  will  be  close  to  the  normal  distri¬ 
bution. 

Let  g  be  the  number  of  secondary  events  NCt)  with  parameters 

A  and  B  let  K  be  a  Poisson  r.v.  with  parameter  y;  let  N  ^  ^  normal 

^  2 
r.v.  with  mean  p  and  variance  o^. 

The  Poisson  character  of  L,  _  for  A  and  B  small  is  described  in 

A»d 

the  following  theorem. 

Theorem  5.1.  For  fixed  a  >  0  ^A.a/A  *  h/2  A  " 

Proof:  The  generating  function  for  is 

o^(v)  •  exp{-(a/A)(l  -  - ' 

As  A  -►  0^  we  find  from  L'Hospital  rule  that  lim  p«,(v)  ■  exp{-»  (1-v)} 

A-^0*  '  ^ 

which  is  the  generating  function  of 

It  is  well  known  that  is  close  to  normal  when  y  is  large. 

Theorem  5.1  then  say  that  the  distribution  will  look  Poisson  when 

2 

E[N(t)]  •  Xpt  /2  is  modest  and  Xt  is  small.  The  asymptotic  normality 
for  the  time  homogeneous  case  will  be  treated  in  the  following  discussion 
for  the  nonhomogeneous  case. 

From  (4.2)  and  (4.11)  we  know  that  N(t)  ■  where  S(t)  > 

Iq  M(T)dT,  with  M(t)  being  the  total  number  of  primary  events  up  to 
time  T  generated  by  a  nonhomogeneous  Poisson  process.  Differentiating 
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(4.13)  with  respect  to  s  at  s  >  0  we  find 

t 

5  *  E[S(t)]  «  /  X(T)(t-T)dT  =  X(t)*t  . 

0 

Differentiating  twice  one  finds  that 


E[S^(t)]  »  X(t)*t^  ♦  {X(t)*t}^. 


(5.1) 


(5.2) 


Hence, 


X(t)*t^ 


(5.3) 


For  (5.1),  (5.2)  and  (5.3),  X(t)  must  be  integrable.  For  the  practical 
cases  of  interest  X(t)  will  be  continuous  and  differentiable.  We 
assume  as  much  smoothness  as  needed  for  what  follows. 

Lemma  5.2.  Let  0.  Then 

d 

(S(t)  -  S^)/o^  as  t 

Proof:  First  note  that 


E[e 


-s(S(t)-§^) 


.2  t 


]  «  exp{^  /  'l'(S,T.t)(t-T)'x(T)dT} 


where 


T(s 


■T,.)  . 


(t-T)V/2 
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Then 


E[e 


-s(S{t)-SJ/o 


‘1 


s^/2  , 

e  '  exp{- 


.2  t 

V/ 

^  0 


[l-VI 


It  is  easy  to  show  that 


1  -  '♦'(-^,T,t)  <  -i  —  ,  for  0  <  s  <  »  . 
°t  ^  °t 

For  s  fixed,  one  then  has 


As  t  the  second  exponent  in  (5.4)  goest  to  lero, 
and  the  normality  follows. 

From  N(t)  ■  *^yS(t)  finds  easily  that 

«  E[N(t)]  -  E(vS(t)]  • 
and 

*  Var[N(t))  • 

We  can  now  prove  the  main  theorm  of  this  section. 

2 

Theorem  5.3.  If  •  as  t-»  ••  then 

(N(t)  -  5  Np  j  . 


:^.T.t)lX(T)(t-T 

t 

- 2 - 

"t 

(5.4) 


(5.5) 

since  t/o^  -»•  0 


(5.6) 


(5.7) 
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Proof:  It  follows  from  (4.12)  that 


and  hence 


-u(N(t)-N.)/t^  -u/C. 

E(e  -  E{exp{-u(l-e  ^)S{t)  ♦  uu§^/C^}] 

-u/C 

-  E[exp{-yCl-e  ') (S(t)-S^) )] 

•E[exp{-uS  (1-e  -  ■^)  ]•  (5.8) 

From  Lemma  S.2,  we  see  that  as  t  -►  »  the  first  expectation  in  (5.8) 

2 

converges  to  u  /2  and  the  second  expectation  to  one.  This  proves  the 
normality. 

Remark  5.4.  The  condition  •  is  not  very  restrictive.  It  holds 

for  a  wide  range  of  rates  X(t).  It  holds,  for  example,  when  X(t)  is 
a  positive  constant  on  a  finite  interval,  or  when  X(t)  ■  t“,  -•  <  o  <  •, 
or  when  X(t)  •  exp  (-yt),  y  >  0.  It  does  not  hold,  for  instance,  for 
X(t)  •  exp  (et),  8  >  0. 

Remark  S.5.  When  the  primary  stream  is  homogeneous  Poisson  we  have 
X(t)  =  X  and 

*  t  , 

5,  -  X  /  (t-t)dt  ■  Xt'/2 

'  0 


and 
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2  ^  ? 

0*  *  ^  /  (t-TD^^dT  -  XtV3  . 

0 

Consequently 

and 

*  u^Xt^/3  ♦  wXt^/2  . 

which  we  already  obtained  in  Section  2.  Clearly, 

3 

>1^.4... 

^t  Xt^/2  ^ 

as  t  -  and  the  asymptotic  normality  is  justified. 
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6.  Practical  Application  of  the  Results 

When  a  new  system  comes  into  use,  its  future  growth  is  often 
uncertain  at  least  to  some  degree.  For  military  systems  growth 
may  depend  on  field  performance  of  the  systems,  on  the  perception 
of  real  costs  and  on  the  uncertainities  in  future  government  funding. 
For  commercial  systems,  market  acceptance,  time  dependent  price  as 
influenced  by  volume  and  competition  contribute  to  the  uncertainty 
of  future  growth.  The  primary  input  rate  \(t)  is  correspondingly 
unknown.  Moreover  the  failure  rate  u  of  a  particular  system  com¬ 
ponent  may  only  emerge  from  experience.  To  avoid  premature  capital 
outlay  costs  and  storage  costs,  scheduled  incremental  acquisition 
and/or  production  will 'often  be  appropriate.  In  a  typical  applica¬ 
tion,  planning  will  be  adaptive  and  ad  hoc. 

The  reader  may  question  the  legitimacy  of  the  assumed  Poisson 
rate  X(t)  for  introduction  of  a  system  into  service.  This  assumption 
is  made  to  reflect  randomness  in  the  arrival  process  and  to  assure 
tractability.  For  a  commercial  system,  such  as  a  commercial  aircraft 
or  large  computer,  where  market  place  acceptance  has  a  strong  role, 
the  arrival  process  has  a  character  similar  to  that  of  rumor  propaga¬ 
tion  Ccf.  Dietz  (1967),  or  Bailey  (1975))  and  the  Poisson  assumption 
seems  indicated.  Even  in  military  contexts  where  a  fixed  schedule 
of  introduction  is  hoped  for,  randomness  in  shipping  times,  training 
times,  etc,  encourages  the  Poisson  assumption.  In  these  few  cases 
where  introduction  at  fixed  intervals  occurs  (cf.  Example  4.S)  the 
formulae  needed  are  obtained  easily. 
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The  adaptive  ad  hoc  planning  will  have  only  one  realization 
(sample  path]  available.  At  time  t  one  will  have  a  history  for 
that  sample  path  of  cumulative  initiations  MCu^.t'),  0  s  t'  <  t, 
to  date.  One  will  also  have  for  each  part  type  A  a  history  of  cumula¬ 
tive  replacements  N^Cu^.t*),  0  s  t'  s  t.  A  residual  number  R^(t)  of 
spare  parts  on  hand  will  be  known.  Suppose  the  failure  rate  g  is 
known.  Let  tj^  ■  t  be  some  future  horizon  time  of  concern.  A 

projected  initiation  rate  X(t')  in  the  period  t  <  t'  <  t^^  is  assumed 
available.  One  can  then  establish  either  from  (3.1)  or  from  the 
result  of  Theorem  5.3  a  probability  of  spare  part  depletion  by  time 
t^  in  the  absence  of  new  spare  parts  production  in  (t,tjj).  The  pro¬ 
duction  quantities  and  the  production  schedule  may  then  be  established 
from  lot  size  techniques  such  as  developed  for  time  varying  deter¬ 
ministic  demand  (cf.  Peterson  and  Silver  (1979,  pp.  300-341)). 

If  the  failure  rate  u  is  not  initially  known,  it  may  be  estimated 
at  time  t  from  cumulative  failure  history,  as  for  instance,  in  Keilson 
and  Sumita  (1982).  In  the  latter  paper  distributed  arrival  times 
corresponding  to  X(t)  are  incorporated. 
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